When do colliding bubbles produce an expanding universe? 
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Abstract 



It is intriguing to consider the possibility that the Big Bang of the stan- 
dard (3+1) dimensional cosmology originated from the collision of two branes 
within a higher dimensional spacetime, leading to the production of a large 
amount of entropy. In this paper we study, subject to certain well-defined 
assumptions, under what conditions such a collision leads to an expanding 
universe. We assume the absence of novel physics, so that ordinary (4+1)- 
dimensional Einstein gravity remains a valid approximation. It is necessary 
that the fifth dimension not become degenerate at the moment of collision. 
First the case of a symmetric collision of infinitely thin branes having a hy- 
perbolic or flat spatial geometry is considered. We find that a symmetric col- 
lision results in a collapsing universe on the final brane unless the pre-existing 
expansion rate in the bulk just prior to the collision is sufficiently large in 
comparison to the momentum transfer in the fifth dimension. Such prior ex- 
pansion may either result from negative spatial curvature or from a positive 
five-dimensional cosmological constant. The relevance of these findings to 
the Colliding Bubble Braneworld Universe scenario is discussed. Finally, re- 
sults from a numerical study of colliding thick- wall branes is presented, which 
confirm the results of the thin-wall approximation. 
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I. INTRODUCTION. 



Recently much work motivated by ideas from string/M-theory has been devoted to the 
study of large extra dimensions [1,2]. Much of this work has been done within the so-called 
braneworld scenarios, where all the low-energy excitations except for gravity are confined to 
a brane, embedded in the simplest case in five- dimensional anti de Sitter space (AdS 5 ) [2]. 
It has been shown that in this type of model, gravity behaves in the large distance limit as 
ordinary (3 + 1) dimensional gravity [2,3]. It has also been shown that in the limit where the 
energy density on the brane is small, the standard four-dimensional cosmological evolution 
on the brane is recovered with calculable corrections that become relevant at high densities 
[4]- 

Given this starting point, there exist a number of ways to proceed toward constructing a 
more complete cosmological model. Perhaps the most conservative route is to postulate an 
epoch of inflationary expansion on the brane. However, when the extra dimensions are large, 
it is not at all clear that such inflation can account for the high degree of smoothness of our 
present universe. A (4+l)-dimensional ("bulk") smoothness problem arises [5,6]. Unless 
some additional mechanism exists for imposing smoothness on the five-dimensional bulk 
spacetime, after inflation on the brane ends, its subsequent evolution into an irregular bulk is 
likely to produce unacceptably large irregularities on the brane. Consequently, considerable 
motivation exists for exploring possibilities for a more complete story for a braneworld 
cosmology, possibly not involving inflation at all. 

A number of models, in particular brane inflation [7], the ekpyrotic universe [8], and 
the colliding bubble braneworld scenario [5] , postulate that the large entropy of the present 
universe has its origin in the collision of branes. In the latter scenario [5,9-11] two expanding 
bubbles collide to form a final brane on which our (3+l)-dimensional FRW universe is 
situated. In this model, the universe starts in a metastable state with a (4+l)-dimensional 
Minkowski or de Sitter geometry. In the former case a natural preferred ground state 
exists. In the latter case, an epoch of old inflation a la Guth rids the spacetime of any 
pre-existing irregularities, effectively establishing an 5*0(5, 1) symmetric initial state with 
calculable fluctuations departing from this symmetry. This initial state decays through 
quantum tunneling by forming bubbles filled with anti de Sitter space [12,13]. If the true 
vacuum is discretely degenerate, the collision of two such bubbles results in the formation 
of a domain wall interpolating between two different such vacua. This is our braneworld. If 
the tension of this local brane is tuned to balance the negative cosmological constant in the 
bulk, the cosmology for an observer situated on this brane is virtually identical to that for 
a Randall-Sundrum model. The decay of one bubble gives a spacetime with its symmetry 
reduced to S0(4, 1). For two colliding bubbles, the symmetry reduces further to S0(3, 1) 
[14,15]. On our braneworld, which has a three-dimensional hyperbolic geometry i?(3), this 
is the symmetry group of spatial manifold at constant cosmic time. If the curvature radius 
today Rq of this manifold is sufficiently large, the universe will appear very nearly flat. 

After the collision, the evolution of the FRW universe is governed by the modified FRW 
equation given in ref. [4]. Given the density on the final brane, which is fixed by stress-energy- 
momentum conservation, the future time evolution of this FRW universe is completely de- 
termined. However, because the Hubble constant on the brane H appears quadratically in 
this equation, it is necessary to determine the sign of H on the final brane. 



2 




fifth dimension 

FIG. 1. Brane Collision Geometry. Two branes collide to form a single final brane sur- 
rounded by pristine bulk spacetime. Here each point represents an manifold of variable 
curvature radius R comprising the three transverse dimensions. 

In this paper, we examine the outcome of symmetric brane collisions, first in the infinitely 
thin brane limit and then numerically for the collision of thick branes. Although this work 
was initially motivated by the Colliding Bubble Braneworld scenario, the calculations and 
techniques developed are more generally applicable. 

The organisation of the paper is as follows. In section II we calculate the outcome of 
collisions of thin branes. We also indicate how these results would be modified when our 
assumptions are relaxed. We show how the same conclusions may be obtained more generally 
by considering the focusing of geodesies as a result of the momentum transfer in the collision. 
In sections III and IV we show how to simulate bubble collisions numerically and present 
the results of our numerical simulations, which agree with the analytic treatment in section 
II. Finally, in section V we present some concluding remarks. 



II. COLLISIONS OF BRANES WITHOUT THICKNESS 

In this section we calculate the outcome of the symmetric collision of two infinitely thin 
branes that join to form a single brane in the final state. It is assumed that the resulting 
brane is surrounded by AdS 5 and that no debris emanates from the collision. Brane collisions 
in the thin- wall limit have been previously considered by a number of authors [5,16-19]. The 
brane collision geometry is sketched in Fig. 1, where only two of the (4+1) dimensions are 
shown. Because of a generalization of Birkhoff 's theorem, the outcome of the collision may 
be analysed by pasting together across the branes the two exterior AdS 5 spacetimes and the 
M 5 spacetime initially between the two colliding branes in such a way that R, the curvature 
radius of the transverse manifold, is continuous across the cuts and that the proper 
time differences on the two sides of the brane agree. Equivalently, it suffices to match R at 
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R = 




R=0 



FIG. 2. Hyperbolic coordinates for describing brane collisions in AdS. Squares such 
as the one illustrated above tile the infinite vertical strip comprising the Penrose diagram for 
maximally extended AdS space. Here four distinct coordinate patches separated by diagonal 
coordinate singularities cover the fundamental region. 



one point and to require that the "Hubble constant" 

B = f T HR\ (i) 

agree on the two sides of the brane. Here r denotes the proper time along the brane. 
H(t) is entirely intrinsic to the brane and fixes the trajectories of the boundaries in the 
bulk spacetimes that have been glued together across the branes. The matching at the 
collision vertex is fixed by the absence of a conical curvature singularity and by stress- 
energy conservation. Note that the equation of state on the brane is necessary to determine 
H; however, for matching at a collision vertex as considered here only knowledge of H is 
required. 

In order to describe the surrounding AdS 5 bulk, having a curvature radius £, it is nec- 
essary to employ coordinates that make manifest the if( 3 ) symmetry of the transverse di- 
mensions. There are two such coordinate systems, covering different patches of AdS 5 and 
separated by coordinate singularities. The type I patches have the line element 

ds 2 = -dt 2 + cos 2 [t/£}d X 2 + I 2 sin 2 [ [t/£]dHf 3) , (2) 

which in the M 5 (£ — > oo) limit becomes 

ds 2 = -dT 2 + dX 2 + T 2 dHf 3) . (3) 

The type II patches have the line element 

ds 2 = -£ 2 sinh 2 [ X /£] dt 2 + d X 2 + f cosh 2 [ X /£] dH 2 {3) , (4) 

which in the M 5 limit becomes the Rindler wedge 
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N 

FIG. 3. Expanding Bubble in Hyperbolic Coordinates. We show the placement of an 
expanding bubble with nucleation center N with respect to the above hyperbolic coordinates. 



We now examine how the conformal diagram of maximally extended AdS 5 is covered 
by such coordinate patches. The conformal diagram consists of an infinite vertical strip 
of finite width tiled by an infinite sequence of square fundamental regions. Fig. 2 shows 
how each of these fundamental regions is covered by type I and type II coordinate patches. 
Note that each point in this figure represents an manifold of variable curvature radius 
R along the three transverse dimensions. The lower horizontal boundary, with R — 0, 
represents an infinite length spacelike geodesic S. The region spanned by the coordinate 
patch of type I just above S is generated by the family of all timelike geodesies emanating 
normally from S. These geodesies are focused in the horizontal direction to join at a point 
situated in the middle of the square, which represents an manifold with R — £. The 
two diagonal boundaries represent coordinate singularities along which this patch joins two 
type II patches. Note that within the type I patches, < R < £. Inside the pair of type II 
patches, £ < R < +oo. Fig. 3 shows how an isolated expanding bubble is situated in relation 
to these coordinates. Since we are interested in bubbles that collide after much expansion 
has taken place, we consider a collision in a type II patch, where R > £. 

A collision may be characterized by two parameters: the curvature radius at collision 
Rcoii, and the boost parameter /3 co a, defined such that the relative velocity of the colliding 
branes is v re i = tanh[2/3 coH ]. Given R cM and j3 C oii, we may match, on the M 5 side with 
geometry (3), onto the trajectory of an expanding bubble of the form, 




(5) 



T = a sinh[/3], 



X = a cosh[/3], 



(6) 



where 



a = 



smh[p coll } 



(7) 
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is the curvature radius of the dS 4 hyperboloid swept out. In terms of the if( 3 ) coordinates, 
the expansion rate of the three transverse dimensions just before the collision is given by 



d , r „ n Id/ , r _\ 1 , r _ cosh|731 
if in = — ln[i?] = -— ln(a sinh[/3]J = - coth[/3] - 



Therefore at the collision point we have 



H in. 



i2 



cosh[/3, 



coii 



For the AdS 5 side of the bulk, we take the line element (4) to obtain the scale factor 

R = £cosh[ x /£}. 



(8) 



(9) 



(10) 



It follows that R co u = £cosh[xcoii/£}- The boost parameter with respect to the AdS 5 rest 
frame PacIS is related to H by 



H = sinh[/3 Ad5 ]^ln(/ cosh[ X /£]) = ^sinh[/3 Ad5 ] tanh[ X /4 



Therefore, just before the collision 



H in = jsmh[PZ s ] tanh[xW4 



It follows that 



cothfe] = 



1 + 



\ m n ) 2 R 2 coU 

The absence of the conical defect at the point of collision implies that 

oout oin o 

PAdS — PAdS ~ Perth 



(11) 



(12) 



(13) 



(14) 



where the hyperbolic angles /3 co ih pAdSi an< l Paos indicate the orientations of the branes with 
respect to the apparent bulk spacetime rest frames, which have been cut and pasted together 
as shown in Fig. 4. It follows that the expansion rate on the final brane just subsequent to 
the collision is given by 



1 i 

H out = -sinh[/^* 5 ] tanh[x co «/^] = - smh[(3% dS - (3 coll \ tanh[x co «/^] 



cosh[/3 coH ] - smh[(3 coll } coth[^ 5 ] 



(15) 



For M 5 initially between the branes, one has 

Hout 



cosh. 2 [P co u} sinh[/3, 



coll\ 



1 + 



£ 2 

^coll 



smh 2 [f3 coU }. 



Rcou £ 

For a bubble of critical radius a, we may re-express H out in the form 



(16) 
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FIG. 4. Hyperbolic angles characterizing the brane collision geometry. In the exterior 
AdS a positive angle corresponds to the direction of an expanding brane. The dotted line in the 
AdS space indicates a contour of constant R, or equivalently the worldline of an observer at rest 
with respect to the AdS bulk. 



H ou t — 




(17) 



We observe that since the term in parenthesis is always negative, the expansion rate for the 
final brane is negative for sufficiently large R co ii- 

We now generalize to the case where the intervening spacetime is dS 5 , having a curvature 
radius i^s, rather than M 5 . In this case, the region initially between the branes is described 
by the metric 



ds 2 = _ dt 2 + cosh 2 [ t /p dS ]d x 2 + f dS smh 2 [t/P dS }dH 2 3) , 
and it follows that 

cosh[/3 coH l _ R: 



Hi, 



v ~dS 



■ cosh[/3 coH ] coth[r/4 5 ] 



R 



coll 



\ 



1 + 



coll 



P 2 

l dS 



(18) 



(19) 



For dS 5 initially between the branes, one has 
cosh 2 [(3 coU } 



H, 



out 



R 



coll 



\ 



rt-coll 

P 2 

l dS 



sinhf/3, 
I 



colli 



\ 



n coll 



sinh 2 [r3 coll ] + cosh 2 \J3 coll ] . (20) 

Pa 



"dS 



For an expanding bubble, R and (3 are not independent. Their interrelation is readily 
described by embedding a bubble trajectory having the geometry of dS A and a curvature 
radius equal to a = Pds sin [a] in dS b of curvature radius Pds (For details see [9]). We construct 
dS 5 through the following embedding in (5+1) dimensional Minkowski space: 
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T = £ dS smh[t/i dS } cosh[£], 

U = £ dS cosh[t/£ dS ] cos[x/4s], 

V = £ dS cosh[t/£ dS ] sm[x/£ dS ], 
X = £ dS smh[t/£ dS } sinh[£] n x , 

Y = £ dS smh[t/£ dS } sinh[£] n y , 

Z = £ dS smh[t/£ dS ] sinh[£] n z , (21) 

having the line element previously given in eqn. (18), where directional cosines. 

dS 4 is embedded as the surface 

U = £ dS cosh[t / £ dS ] cos[x/£d,s] = £dscos[a] = (constant). (22) 

It follows that the velocity with respect to the rest frame of the metric in eqn. (18) is 

4. uia i uu//> i d * uuie i tanh [*/4s] R coU cos[a] 

v = tanh [(3 coU \ = cosh [t £ dS \ — = cosh [t/£ dS - — L , J = . 23 

dt tan[x/4s] yJR^ + a 2 

where R = £ d s smh[t/£ d s] so that 1 

C0Sh[M= ^^ = , sinh[fall= g^sa . ( 24 ) 



It follows that 



-dS 



As R coU -> oo, 



Consequently, when 



or equivalently when 



^ sin 2 [a] £ dS sm[a]U 2 a 2 ' W 



/ a 2 \ 1/2 

cos[a] 1 + - <1, (27) 



;-£)H)<^ 



1 Note that as i? co /z — >• oo, cosh[/3 co ;;] — >• csc[<j], which is finite rather than infinite. Because of the 
inflation of the intervening spacetime, colliding bubbles expanding in dS 5 have a finite terminal 
velocity. 
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a situation occurs where the final brane is always expanding, for arbitrary large radius of 
curvature R co ii of the surface of collision. 

The condition in eqn. (28) is equivalent to the condition 

c?>e dS -e. (29) 

Finally, we remark on what happens when some of our simplifying assumptions made 
to facilitate this calculation are relaxed. In our calculation we assumed that all the energy 
available in the collision is deposited on the final brane. One may ask what happens when 
some of this energy is carried away by other branes or by radiation. We may also ask what 
happens when the branes are not infinitely thin. We retain the assumption of a symmetric 
collision. 

First, it is possible to show that if energy, or more precisely momentum, escapes from the 
collision away from the central final brane, H out becomes smaller or more negative. This can 
be seen by considering the family of geodesies in the midplane of the collision, which after 
the collision become the fiducial observers on the FRW universe on the final brane. These 
fiducial observers can be continued back into the past prior to the collision into the initial 
intervening bulk spacetime, where they have an initial Hubble constant H be f orej cen t e r, which 
suffers a discrete jump at the collision. This jump is proportional to the impulse in T| at 
the collision, which in turn is proportional to the momentum transfer across the collision. 
More explicitly, the Einstein equations in five dimensions give on the 4-D plane of symmetry 

R R 2 « _ 87tG( 5 ) 

R R 2 R 2 3 5 

Here G( 5 ) is the five dimensional Newton constant, k/R 2 the spatial curvature term, and T| 
the "pressure" along the fifth dimension. Assuming a sharp impulse, we obtain 

(tt tt \ {Rafter ~ Rbefore) 8ttG (5) f+ e 5 87tG( 5 ) 5 

[H after - H before ) cent£r = = — dr T 5 = — J 5 , (31) 

where e is taken to be infinitesimal. If follows from stress-energy conservation that 

^5 = Pafter ~ Pbefore' (32) 

where p 5 is taken positive in the sense away from the collision. Eqn. (32) may be verified 
by surrounding the collision by an infinitesimal pillbox and integrating T^ 5 over its surface. 

The integral If is therefore minimized in the calculation presented above, where the 
collision resembles an 'inelastic' collision with nothing bouncing off. If the branes do not 
completely stick to each other, Paf ter is increased, consequently decreasing H a f ter . Note that 
for thick branes the same analysis may be repeated in the midplane of the brane. 

III. NUMERICAL CALCULATIONS. 

We now present some numerical simulations of collisions of expanding bubbles to explore 
what occurs in the thick-wall case. To study the gravitational effects of the bubble collision 
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(30) 



V[4>] 




4> 



AdS\ AdS2 



FIG. 5. Scalar Field Potential. The initial state for the scalar field is a metastable minimum 
at 4> = where the spacetime geometry is either M 5 or dS 5 . The field decays to either of the two 
AdS true vacua. 

we will concentrate on a model with a single scalar field coupled to gravity described by the 
action, 2 



where the self interacting potential for the scalar field is of the type drawn in Fig. 5. This 
simple potential has the correct properties to describe the colliding bubble braneworld sce- 
nario. We assume an initial state where the scalar field is trapped in a metastable state 
at (p = 0, where cosmological constant is either positive or zero so that the bulk spacetime 
outside the bubbles is either de Sitter or Minkowski space, respectively. The field tunnels 
to the true vacuum by the nucleation of bubbles filled with anti de Sitter space. We first 
discuss the situation with only one bubble, so that the symmetry of the problem is 5*0(4, 1) 
and the equations for the instanton describing the decay of the false vacuum are 




(33) 




(34) 



and 




(35) 



where the Euclidean metric of the instanton is 



2 We set 8vrG, 



1 unless otherwise stated. 
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Region II 



FIG. 6. Spacetime diagram of the bubble nucleation. Region I, in the interior of the 
bubble is an open FRW universe. Region II is the exterior of the bubble. The scalar field <j> is 
constant along the lines shown in this picture. 



dr E = da 2 + b 2 (a) df E + cos 2 [t E ]dQ 2 3) 



(36) 



and the prime denotes differentiation with respect to a. Upon analytic continuation of t E 
to t = it e, the line element becomes 



df L 



da 2 + b 2 (a) -dP + cosh 2 [t]dQ 



(3) 



(37) 



describing the instanton solution in region II of Fig. 6. 

This solution may be continued into region I (i.e., the interior of the forward lightcone 
of the nucleation center) by solving the coupled equations 



and 



+ 4 



-V, 



1 • 



2 + V[<t>] 



+ 



with line element 



ds 2 = -dr 2 + a 2 (r) [d X 2 + sinh 2 [ X }dVl 



(38) 



(39) 



(40) 



This system of equations describes the interior of the bubble as an open (hyperbolic) 
universe filled with a scalar field <fi and evolving with a scale factor given by a(r) [13]. Due 
to the negative value of the cosmological constant inside of the bubble, the initial expansion 
of this open universe slows down and eventually reverses, leading to a collapsing phase. In 
the thin- wall limit, = (j) tv exactly in region I and the geometry is precisely that of perfect 
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anti de Sitter space, with a(r) = H^ 1 sm[Hr]. In this idealized case, r = nH^ 1 corresponds 
to a coordinate singularity. However, with any non-vanishing bubble wall thickness, the 
wall always has a decaying tail, so that at the lightcone <f> is always slightly displaced from 
its true vacuum value. For the case of an expanding interior geometry, the deviation from 
(ptv oscillates and decays, with becoming increasingly close to (f) tv as one passes farther 
into the bubble interior. However, when the expansion slows down eventually turning into 
a collapse, the deviation from tp t v has two modes: one regular and the other divergent as 
the "Big Crunch" of the collapse is approached. The "Big Crunch" is a backward directed 
lightcone when the gravitational backreaction of these modes is ignored. Unless the potential 
is finely tuned with infinite precision, at least some amplitude of the divergent mode is 
always present. This divergent mode leads to a spacelike coordinate singularity [13,20], akin 
to the initial singularity of a homogeneous and isotropic universe with a hyperbolic spatial 
geometry. 

The singularity forms as the result of the extremely finely tuned focusing of the scalar field 
toward the same point in a large region of spacetime. The role of the negative cosmological 
constant is subsidiary. It is important only during the middle phase of the evolution between 
the two lightcones, where it acts as a sort of perfect lens without any aberration at all. Once 
the initially divergent field mode has been refocused, the negative cosmological constant is 
no longer needed. The formation of the singularity would occur, even in Minkowski space, 
for example, where all but the last moments of the collapse are described by a sort of upside 
down Milne universe. In the presence of sufficiently large quantum fluctuations exceeding a 
certain threshold, away from the semi-classical (h — > 0) limit, the degraded quality of the 
focusing presumably allows the formation of a singularity to be avoided. 

The original motivation for carrying out these simulations was to study this singularity 
in the presence of a collision with another bubble and its consequences for the colliding 
braneworld scenario (See [11]). The argument was that the collision of scalar field bubbles 
would most likely produce some radiation which would propagate in the background space- 
time described above. This radiation could however spoil the focusing of a certain part of the 
singularity, thus opening a window into infinite AdS space. This is what we set out to study 
by doing numerical simulations of the colliding bubbles in the scalar field model described 
above. The simulations performed, however, indicate the presence of another singularity 
forming in the future of the collision region. This implies that the 4d universe within the 
plane of symmetry that would be our braneworld would rapidly collapse. These numerical 
results led us to investigate this issue further and this paper reports on this line of work. 

A. Fixed Minkowski background 

In order to gain a better understanding of what to expect from the simulation of the 
bubble collision we first simulated a collision with the background fixed to Minkowski space. 
We first find the solution for the scalar field from the instanton equation (34), fixing the 
background metric to be Minkowski, in other words with b(a) = a. We then use this solution 
to generate the initial state with two bubbles, by just pasting together two of these solutions 
at some distance d and reversing the sign of <fi for one of them. We expect this initial data to 
be a good approximation to the real situation of the nucleation of two bubbles at the same 
time at sufficiently large separation d. To evolve this forward in time, we make use of the 
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symmetry of the problem by choosing the hyperbolic slicing of Minkowski space given by 



ds 2 = -dt 2 + dx 2 + t 2 dHf 3) , (41) 

where dH 2 ^ = d£ 2 + sinh[^]dQ 2 2 y This is just Minkowski space expressed as the direct 
product of a Milne universe in one dimension less with the line — oo < x < +00 of a spatial 
character. It is along this direction x where the two bubbles are separated by a distance d 
initially. In this coordinate system the equation of motion for the scalar field becomes 

d 2 <P 3d<P d 2 <P_ 

W + lTt-^-~ v ^ (42) 

Note that even though this coordinate system does not cover the whole spacetime, it allows 
us to study the region of the collision between the two bubbles. 

In Fig. 7 we plot the isosurfaces of constant scalar field in the (x, t) plane. Each point in 
this 2D picture represents a hyperboloid of radius t on which the scalar field is constant. 
We see from this simulation that the bubbles start at rest with an initial radius a and 
rapidly accelerate toward one another until finally colliding with a considerable amount of 
kinetic energy at the plane of symmetry x — 0. At the collision the two bubbles bounce off 
one another while continuing to expand but the attractive force between them bring them 
together again and again until they finally settle down to a stable configuration of a kind of 
"thick" domain wall or brane. This is the braneworld where we live. 

One important point here is boundary conditions. We have imposed symmetric boundary 
conditions for the field <fi and this allows for the radiation to be reflected on the boundaries 
and enter the region of interest again. This is not desirable and we will see how we do not 
have to impose any boundary conditions of this type when we include gravity. 



B. Including gravity 

We now study the same problem including gravity. As before, it is most convenient 
to employ a system of coordinates that maximally exploits the 5*0(3, 1) symmetry of the 
problem. To cover all the spacetime in this way would require several coordinate patches, 
separated by boundaries across which the geometry of the orbit of a point under the action 
of SO(3, 1) changes from to dS( 3 y However, for our purposes, because our primary 
interest is the region of collision and the brane between the bubbles arising from this region, 
it suffices to consider one such patch, described by a metric of the form 

ds 2 = -f 2 (u, v) du dv + r 2 (u, v) dHf 3) . (43) 

To further fix the gauge, let us additionally restrict the coordinate patch to the region 
(u + v) > and require that r = on the line (u + v) = 0. The gauge may be completely 
fixed by requiring that f 2 (u, —u) — 1. In this gauge, we may rename cr+ = u, <r_ = —v, so 
that the metric becomes 

ds 2 = +/ 2 (<7-, 0+) da + da_ + r 2 ((j_, 0+) dH 2 3) (44) 

where, physically, the coordinates a + and o"_ of a point P indicate the intersections of the 
right-moving and left-moving null geodesies passing through P with the r = line, where a 
represents the physical distance along this line. 
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FIG. 7. Isosurfaces of constant scalar field <fi for a collision of two bubbles in the (x, t) 
plane of the 5D-Minkowski space. The bubbles collide and bounce off each other emitting 
radiation at each encounter till finally joining to form a thick brane that interpolates between the 
two AdS minima of the scalar potential. 



We now show how to describe the general expanding bubble solution, endowed with an 
5*0(4, 1) symmetry in terms of the SO(3, 1) symmetric coordinates just described. In terms 
of the natural coordinates of the larger SO (A, 1) symmetry, the solution is described using 
two coordinate patches, a region I (inside the forward lightcone of the nucleation center), 
with the metric 



ds 2 = -dr 2 + a 2 (r) dH 2 (A) 

= -dr 2 + a 2 (r) d^ + sinh 2 [x]^ 2 3) 
= -dr 2 + a 2 (r) [d£ 2 + cosh 2 [£]dff ( 2 3) ] , 

and a region II, outside the lightcone, with the metric 

ds 2 = da 2 + b 2 {a) d(dS) 2 {A) 



(45) 



= da 2 + b 2 (a) 
= da 2 + b 2 (a) 



-dP + cosh 2 [t]dQ 2 



(3) 



-di 2 + smh 2 [i\dHf 3) 



(46) 



The coordinates on the last line of the previous equations cover only a subregion of the region 
covered by the coordinates of the preceding line; however, for studying bubble collisions this 
does not pose a problem. 

To treat both regions within a common framework, it is useful to transform to null 
coordinates, denning the transformed variables 





i r 


[ T dr 


a(r) t 


Jo 



ln[r] +g(r), 



(47) 
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and 



ln[o] + f 
Jo 



do 



6(a) 



Ha} + f(a), 



(48) 



so that eqns. (45) and (46) become 



ds 2 = a 2 (r) [-d(t + £)d(t - + cosh 2 [£]dff ( 2 3) 
= a 2 (t at) ) + cosh 2 [£]dif 2 3) 



and 



= a 2 


(4*= 




—dt + dt_ 


+ cosh 2 


■t+ 


d# ( 2 3) 




) 


2 


ds 2 = b 2 ( 


er) d(s 4- 


t )d(s - i ) 


+ sinh 2 [t]d# 2 3) 




= b 2 


/s+ + s_ 
V 2 


■) 


—ds + ds_ 


. + sinh 2 


"s+ - s_" 
2 


dff ( 2 3) 



(49) 



(50) 



respectively. To patch together the two regions, it is convenient to rescale the null variables 
back to r and o, so that the metric is 



ds 2 = a 2 {r av ) 



—dr_dr^ 
a(T + )a(r- 



+ cosh 2 



t+-t. 



and 



ds =b (a av ) 



do -do. 



+ 



[b(o + )b(o-) 



+ sinh 



s + — s_ 



dH (3) 



dH (3) 



(51) 



(52) 



Note that the surface of the lightcone corresponds to the limit r + =(fixed), t_ — > 0+, and 
o + = (fixed), <r_ — > 0+ in the respective regions inside and outside the lightcone. Conse- 
quently, the lightcone is parameterized equally well by both o + and r + . The correspondence 
between these variables is established by matching the curvature radius of the H^) dimen- 
sions by setting 



lim a 2 (T a v) cosh 
From eqn. (47), it follows that 



U-t- 



lim b 2 (o av ) sinh 2 



r av e 9( - Tav) = e tav = e {t++t - )/2 = v /7vr7e [9(r+)+9(T - )1/2 . 
However, as t_ — > 0, <?(t_) — > 0, and 

Too = JwFZePW 2 , 



(53) 



(54) 



(55) 



and since a(r) = r + 0(t 3 ) (for else the curvature on the lightcone would be singular), it 
follows that 



a(r av ) = 0v^e^+)/ 2 . 



(56) 
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Similarly, in the same limit 



cosh 



■t+ - t. 



cxp 



(t+-t_) « — e 9(r+) . 



(57) 



Hence 



lim a 2 (r a „) cosh 2 

T_^0+ 



T 2 e 2 5 (r + )_ 



Likewise, the limit on the right-hand side is equal to a^_e 2 ^ a+ \ It follows that 



r+e 9(r + ) = (J+e /(-+) ) 



(58) 



(59) 



or 



(60) 



By matching in this way, it is possible to include both the bubble interior and exterior 
within the same coordinate patch. a+ and cr_ are allowed to take both positive and negative 
values, subject to the restriction er_ < a + . The range cr_ < < er + covers the interior of the 
lightcone. Geometrically, the transformations (59) and (60) map the a line into the r axis 
such that null geodesies link the two axes. 

If we write the metric in the form 



dr 



f 2 (a_, a + )da_da + + r 2 (o"_, a + )dHf 3) , 



(61) 



it follows that 



b 2 (cr av ) 
b(a+)b(a-)i 

a (jay) 

6(a+W-^-)' 

b 2 ((Tav) 

b(-a+)b(-a-)i 



< a- < <t+, 
o- < < a + , 
cr_ < (7+ < 0, 



(62) 



and 



rV_,cr + ) = < 



b 2 (a av ) sinh 2 
a 2 (T av ) cosh 2 
b 2 (a av ) sinh 2 



2 

-s(cr + )-s(-a-) 
2 



< a- < cr+, 
o- < < cr+, 
o_ < a . < 0. 



(63) 



By considering the appropriate limits, it may be verified that both / 2 and r 2 are smooth 
across the lightcone. 

We are now ready to rewrite the whole spacetime solution of the bubble nucleation in 
terms of the null coordinate system natural to this problem. This will be crucial in the next 
section in order to obtain the initial data for the numerical simulation of the bubble collision 
in terms of the one bubble solutions. 
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100 200 300 400 500 

FIG. 8. Approximate region of the spacetime simulated. We show the scalar field 

solution from the flat background case for reference. It is this region within the black square that 

we plot in Figs. 10 and 11. 

1. Equations of motion for a doubly null coordinate system. 

In the doubly null coordinate system given by the metric 

df = -f 2 (u, v) du dv + r 2 (u, v) dHf 3) , (64) 
we can write Einstein's equations for a minimally coupled scalar field in the following way 

r uu f ~ 2f u r u + %-fct>l = 0, (65) 
r w f - 2f v r v + ^rf4> 2 v = 0, (66) 
2rr uv + Ar u r v -f- V fV = 0, (67) 

2r-V™ + f-'f uv + r-\ u r v - f~ 2 f u f v - \fr~ 2 + (^L -Vj=0, (68) 

and the scalar field equation becomes 

3 1 

(puv + -r~\r u <p v + r v (f) u ) + -fVj = 0. (69) 

In order to integrate these equations of motion numerically, we followed the techniques 
developed by Hamade and Stewart [21] with appropriate modifications. We transformed the 
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FIG. 9. Spacetime diagram for the bubble collision. The initial data for the simulation 
is placed on two null surfaces Ml and Mr. 

problem into the system of first order differential equations described in the Appendix. We 
find it convenient to specify initial data on a pair of null surfaces Ml and Mr at right angles 
with respect to each other, chosen so that the bubble collision occurs in the center of the 
resulting future diamond. See Fig. 8. 

In our case we use the prescription given in the previous section to obtain the form of the 
functions f(u, v), r(u, v) and (f>(u, v) and their derivatives on the null surfaces Ml,r in Fig. 9 
from the known functions of the metric a(r) and b(a) as well as from the scalar field in the 
one bubble case. We follow the same algorithm as in the fixed Minkowski space case and 
copy the same single bubble data on both branches of M. It is clear from our choice of this 
surface that this is a good approximation if we are mainly interested in what happens around 
the collision region, since the existence of the second bubble would not greatly change the 
data on this surface. 

IV. NUMERICAL CHECKS AND RESULTS 

To verify the correctness of our code, we first simulated the evolution of a single bubble 
to enable comparison to the solution obtained from eqns. (34)-(35) and (38)-(40). This is a 
nontrivial check since it tests all the technically difficult steps required to go from the two 
effectively one-dimensional solutions inside and outside of the bubble 3 to its description in 
terms of the new variables in the (u, v) plane. This check shows that the code based on 
eqns. (65)-(69) reproduces the known solutions for the single bubble case (even close to the 
singularity inside of the bubble) more and more accurately as we increase the fineness of our 
grid. 



3 We can obtain the solutions for region I and II separately by solving the one-dimensional 
eqns. (34)-(35) and (38)-(40). This is due to the extra symmetry present in the one-bubble 
case. 
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FIG. 10. Isosurfaces of constant r(u, v ) for the collision of two bubbles. 

A. Numerical results 

We have simulated the collision of two bubbles with several potentials and different values 
of the grid lattice step and the results are qualitatively the same. The main conclusion is 
illustrated in Fig. 10. Here we show contour plots of the function r(u,v) indicating the 
size of the transverse dimensions. We see that r increases linearly along the symmetry axis 
u = v in accord with the Milne-type description of Minkowski space explained above. This 
continues to hold until the moment of the collision, after which r rapidly declines finally 
vanishing at a singularity. If we just imagine the freely falling observer following along a 
u = v geodesic, she experiences a rapid collapse of the three dimensional transverse space 
dimensions finally crashing into a Big Crunch singularity. Fig. 10 shows that this singularity 
extends in a spacelike manner away from the would-be brane. 

We plot in Fig. 11 the isosurfaces of constant 0. The singularity appears here as a 
divergence in the scalar field along the spacelike surface S. 

These numerical results agree with the analytical results presented earlier. The source 
for the positive pressure in this case is of course the region around the collision where the 
gradient for the scalar field grows from a very small value at the beginning of the simulation. 
We have simulated only the collision of bubbles expanding in Minkowski space. 

V. DISCUSSION 

In this paper we have studied the collisions of (3 + l)-dimensional branes embedded in 
a (4 + l)-dimensional bulk spacetime from two points of view: First, we investigated the 
collisions of thin shells analytically. Then we studied numerically the collisions of thick 
bubbles with gravity included using a two-dimensional doubly null lattice simulation. We 
found that in many cases the outcome of the collision is gravitational collapse to a "Big 
Crunch" singularity on the final brane, although some situations where this collapse is 
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averted were found. The thick-wall numerical simulations suggested that the singularity is 
not confined to the three transverse dimensions, but rather extends in a spacelike manner 
in the fifth dimension. 

In a symmetric collision the tendency toward collapse can be understood as follows. By 
continuing the fiducial observers of the final brane into the past, we can formulate a "jump 
condition" at the bounce, as given in eqn (31). Since for colliding positive energy branes 
the impulse If is always positive, expansion of the final brane (H a f ter > 0) requires a pre- 
existing expansion in the bulk (i.e. Hb e f ore < 0). This impulse always acts in the sense of 
slowing down, if not reversing, the pre-existing bulk expansion. In some cases, Hbefore is 
large enough so that H a f ter is positive. 

Our results for the collision of the thin-wall branes may be demonstrated graphically 
as in Figs. 12 and 13. In Fig. 12 we show that two spatially flat colliding branes initially 
separated by M 5 always produce a collapsing final brane. The dashed lines indicate contours 
of constant size R for the transverse dimensions. This scale is such that R for neighboring 
contours differs by a constant factor. In this case, because R is constant in the M 5 initially 
between the branes, these contours must be parallel to the incident branes. As seen in the 
figure, R must contract in the forward time direction along the final brane. In Fig. 13 we 
illustrate how pre-existing expansion in the bulk can alter this conclusion. Here we show a 
final brane that is expanding. 

We conclude with the following comments: 

1. Our analysis supposes the absence of any new physics beyond conventional classical 
five-dimensional Einstein gravity. It is possible that new physics could alter the conclusions 
of our analysis. For example, the absence of a deficit angle at the collision, as expressed 
in eqn. (14), can be shown to hold for a singular distribution of matter of co-dimension 
one in conventional gravity. New physics or higher order terms, however, could alter this 
conclusion. In the ekpyrotic models where the fifth dimension becomes degenerate at the 
instant of collision, new physics is required and our analysis cannot be applied. 



20 



AdS 5 : 



' contraction 



contraction 




FIG. 12. Collision of two positive energy M 4 branes initially separated by M 5 . The 

diagram above indicates a visual representation of why flat branes with five-dimensional Minkowski 
space initially intervening always collide to form a collapsing final brane. The dashed lines indicate 
contours of constant R, the R of neighboring contours differing by a fixed ratio. Since the branes are 
of positive tension, R decreases away from the incident branes. Because R is uniform throughout 
the M 5 space at the bottom, these contours must be parallel to the incident branes. Consequently, 
R on the final brane must decrease in the forward time direction, indicating collapse. 



2. When the bulk spacetime initially between the colliding branes is M 5 , as indicated in 
eqn. (17), whenever R/a is large there is always collapse on the final brane. For the Colliding 
Bubble Scenario, small values of R/a are not acceptable because fl C oii would be significantly 
smaller than unity, leading to an empty universe today. 

3. For dS 5 initially between the branes, if a 2 > l 2 dS — I 2 there is always expansion for 
arbitrary large R co u- Consequently, for Ids ~ I the final brane will always be expanding. 
However, since a < las, Rcoii/lds would have to be enormous. In reference [22] it was shown 
that no matter how small the nucleation rate is, collisions between large bubbles (with 
Rcoii 3> Ids) would have a strong power law suppression and would therefore be extremely 
unlikely. 
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FIG. 13. Collision resulting in expansion on the final brane. The conclusion displayed 
in the previous figure that the final brane is alway collapsing may be circumvented when there 
is expansion in the bulk prior to the collision, as shown above. This expansion may occur either 
because of a hyperbolic geometry, where M 5 has a Milne geometry, or as the result of a positive 
bulk cosmological constant so that dS 5 takes the place of M 5 . 
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VII. APPENDIX 

A. Equations of motion for the numerical code. 

This Appendix closely follows ref. [21]. We implement a first order algorithm to integrate 
numerically the system of equations given in Section III for the bubble collision. Using the 



following definitions 

•=i/l* v -k v - v -f- (70) 

q = s v , p = s u , j = r u , g = r v , (71) 

c= f j, d=^, (72) 

fi=jq + gp, \ = j g -^f 2 , (73) 
we can rewrite the second order equations of Section III as 

rq u + \v + rfX = 0, rp v + ^ + rfV c = 0, (74) 

rjv + 2A - 2r 2 f 2 V c = 0, rg u + 2A - 2r 2 f 2 V c = 0, (75) 

r 2 c v - 3A + \qr 2 + r 2 f 2 V c = 0, r 2 d u - 3A + \qr 2 + r 2 f 2 V c = 0, (76) 

j u - 2cj + rp 2 = 0, (77) 

g v -2dg + rq 2 = 0. (78) 



These equations must be supplemented, following the prescription given in the main text, 
with initial data on the null surfaces Ml and Mr. 
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